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Abstract 

Motivated by recent interest in permutation arrays, we introduce and in- 
vestigate the more general concept of frequency permutation arrays (FPAs). 
An FPA of length n = mX and distance d is a set T of multipermutations on 
a multiset of m symbols, each repeated with frequency A, such that the Ham- 
ming distance between any distinct x,y G T is at least d. Such arrays have 
potential applications in powerline communication. In this paper, we establish 
basic properties of FPAs, and provide direct constructions for FPAs using a 
range of combinatorial objects, including polynomials over finite fields, com- 
binatorial designs, and codes. We also provide recursive constructions, and 
give bounds for the maximum size of such arrays. 

1 Introduction 

As indicated in j4j and |7|, permutation arrays arise in the study of permutation 
codes, which in turn have a natural applicability to powerhne communications. An 
electric power hne may be used to transmit information in addition to electric power, 
by modulating its frequency to form a set of n close frequencies. These small vari- 
ations may be decoded as symbols at the receiver. Steps must be taken to ensure 
that this information transmission does not interfere with the hue's primary func- 
tion of power transmission, and for this reason block coding is used (codewords of 
fixed length). A code is a constant composition code if each codeword, of length n, 
has precisely occurrences of the i-th symbol, where the are positive integers 
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satisfying ^ = n. (Here, the i-th symbol corresponds to the i-th frequency.) 
Various tradeoffs must be made between the competing goals of addressing noise 
problems and the requirement of a constant power envelope. 

One approach is to choose ri = r2 = • • • = r„ = 1, in which case each codeword 
is a permutation on n symbols. An (ri, d) permutation array, usually denoted by 
PA{n, (i), is a set of permutations of n symbols with the property that the Hamming 
distance between any two distinct permutations in the set is at least d. Permutation 
arrays are important not only in powerline communications as described above; they 
have also been applied in the design of block ciphers; see P . 

In this paper, we introduce a generalization of permutation arrays, which we 
call frequency permutation arrays. These arise from the constant composition codes 
when we take ri = r2 = ■ ■ ■ = = A, for some A such that n = mX. When A = 1, 
this reduces to the permutation case studied in 4^ and . We present various results 
and constructions for frequency permutation arrays, many of which have well-known 
permutation array results as special cases. There is a strong connection with recent 
work on constant composition codes such as [H] and [TT] . 

2 Frequency permutation arrays 

We consider rearrangements of the n-element set {1, . . . , 1, 2, . . . , 2, . . . , m, . . . , m}, 
where n = mX {n, m, A G N) and each of the m distinct symbols occurs exactly 
A times. When A = 1, the set of all such permutations is the symmetric group 
Sn of permutations on n symbols. In the general case, these rearrangements are 
multipermutations on the multiset {1, . . . , 1, 2, . . . , 2, . . . , m, . . . , m} (each symbol 
occuring A times); we shall call them X-permutations. 

Definition 2.1 Two distinct X-permutations cr = si . . . s„, r = ti . . . t„ have dis- 
tance d{a, t) = d if they disagree in d entries, i.e. if \{i '■ Si ^ ti}\ = d. 
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This is the Hamming distance famihar from coding theory. In the case when 
A = 1, two permutations a,T E Sn have distance d if crr~^ has exactly n — d fixed 
points. 

Definition 2.2 A permutation array of length n and minimum distance d, denoted 
by PA{n, d), is a subset T of Sn such that the distance between any two members of 
T is at least d. A PA{n, d) may be viewed as an s x n array whose rows are the s 
permutations ofT in image form; taken pairwise, any two distinct rows differ in at 
least d positions. The maximum possible size of a PA{n, d) is denoted by M{n, d). 

We define a frequency permutation array as follows. 

Definition 2.3 A frequency permutation array of length n = mX, frequency A and 
minimum distance d, denoted by FPAx{n,d), is a set T of X-permutations (multi- 
permutations of the multiset {1, . . . , 1, 2 . . . , 2, . . . , m, . . . , m}, each symbol repeated 
X times ), with the property that the distance between any two members of T is at 
least d. Equivalently, an FPAx{n,d) is an s x n array whose s rows consist of m 
distinct symbols, each repeated exactly X times, such that the distance between any 
two rows is at least d. 

Thus an FPAi{n, d) is simply a PA{n, d). We let Mx{n, d) denote the maximum 
possible number of rows that can exist in any FPAx{n, d); then Mi{n, d) = M{n, d). 

Example 2.4 An FPA3{6,A) of size 4 is given by 

1110 
10 110 
^"010101 
10 11 

We first establish basic properties of frequency permutation arrays. Various basic 
results on permutation arrays (for example from ^ and [7]) appear as special cases 
of these results. 

3 



Theorem 2.5 Let n = mX. Then 
(%) Mx{n,2) = ^; 

(a) M\{n, n) = m; 

(zzz) Mx{n, d) > Mxin - X,d), Mx{n, d + 1); 
(iv) If Til = mXi and n2 = mX2, then 

Mx^+\^{ni + n2,di + c/2) > mm{Mx^{ni,di), Mx^{n2,d2)}. 
In particular, M2x{'2n,2d) > Mx{n,d). 



Moreover, for any divisor I of X, Mx{n,d) < jMi(n,d). 

Proof (i) Since two distinct multipermutations must differ in at least two entries, 
Mx{n, 2) is the number of distinct A-permutations. There are (^) . ("^'^) • • • • (^"^("^^^^-^^ 

choices for each A-permutation, i.e. {"^'x) ■{'^"^'"x^'^) ■ ■ ■ - (a) ~ {^Tp' ^^^^ multipermu- 
tations in total. 

(ii) Since there are at most m choices for the symbol in the first position of a 
A-permutation in an FPAx{n, n), we have Mx{n, n) < m. Take m blocks comprising 
A copies of each symbol: {0 ... 0}, {1, . . . , 1}, . . . , {m — 1, . . . , m — 1}; applying an 
m-cycle to these blocks yields m A-permutations, all of pairwise distance n. 

(iii) Adding A copies of some new symbol to each row of an FPAx{n — A, c?) 
yields an FPAx{n,d); the second observation is immediate from the definition. 

(iv) Juxtaposing an FPAxj^{ni,di) and an FPAa2(^2, (^2) yields an FPAx{n,d) 
with A = Ai + A2, "n, = ni -|- ?T,2 and d = di + ^2- 

(v) This is proved in Theorem 14.11 from the size of a PA[n,d) is bounded 



(v) Mx{n,d)<^;Mx{n,d)< 




above by 




□ 
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For any A-permutation a, the sphere with centre a and radius r is defined to be 
the set of all A-permutations with distance at most r from a. We denote its volume 
by Vx{n,r). 

Lemma 2.6 Let n = mX. Then 

^ ^^tE;^^^© ... Q(-i)'f e-^{n..(^)}-. 

where Lk{x) is the kth Laguerre polynomial. Here the inner sum runs over P{k) = 
{{ki, . . . ,kt]ri, . . . ,rs)}, the set of all partitions ki + ■ ■ ■ + kt of k ^ N into positive 
integers 1 < fcj < A, where the set {ki, . . . , kf} consists ofrj occurrences of value ki. 
(j = 1, . . . , s), with 1 < ki^ < A, 1 < < t and ri + • ■ ■ + = t. 

Proof Let a be any A-permutation of length n. The set of A-permutations at 
distance k from a is obtained by taking each fc-entry subset of a, and deranging its 
entries. By a result obtained in [12j, and reproved in |3], the number of derangements 
of a sequence composed of ni objects of type 1, n2 objects of type 2,. . ., Uf objects 
of type t (i.e. permutations in which no object occupies a site originally occupied 
by an object of the same type) is given by 

poo J' 

D{n,,...,nt) = {-lf / e-^{r[L„^,(x)}rfx, 

where ni + ■ ■ ■ + nt = N . The result follows upon applying this theorem to each 
fc-element subset of a. For any A-permutation a, and any partition ki + ■ ■ ■ + kt 
of G N into positive integers l<A;j<A(l<t< m), we count the number 
of fc-subsets comprising ki occurrences of symbol si, k2 occurrences of symbol S2, 
. . ., kt occurrences of symbol St- Suppose the set {fci, . . . , fcj} consists of ri occur- 
rences of value fcj^, . . ., Ts occurrences of value fcj^, where ri + ■ ■ ■ + Ts = t. There 
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are (™) {"^^p) ■■■('" ^'=i = ^7(m-t)! choices for symbols Si, . . . , Sj. For each 
choice, there are (^J (^^) . . . (^J subsets of a in which elements occur with appro- 
priate frequency. □ 
A covering argument yields the following lower bound for Mx{n, d), an analogue 
of the Gilbert- Varshamov bound in coding theory, while a sphere-packing argument 
yields an upper bound, analogous to the Hamming bound for coding. 

Theorem 2.7 

We remark in passing that a useful upper bound for the maximum size of general 
constant-composition codes (CCCs) has recently been presented in [12] and has been 
further developed in jTTj. However, for a CCC in which all symbols of a codeword 
occur with equal frequency A (the situation corresponding to FPAs), this upper 
bound essentially reduces to the Plotkin bound, Mx{n,d) < , which is valid 

only when d > n ~ X. Since the direct constructions presented in this paper have 
minimum distance less than or equal to n — A, the bound is of limited applicability 
in our setting. 

3 Direct constructions 

It is known that permutation arrays may be constructed using latin squares (see 
|lj and |ini)- Frequency permutation arrays are related to frequency squares as 
permutation arrays are to latin squares, and this connection may be exploited to 
obtain a construction for FPAs. 

Recall that a latin square of order n is an n x n array in which n distinct 
symbols are arranged so that each symbol occurs once in each row and column. 
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Two latin squares Li and L2 of the same order n are said to be orthogonal if, when 
superimposed, each of the possible ordered pairs occurs exactly once. A set 
{Li,L2, . . . ,Lt} of t > 2 latin squares is said to be mutually orthogonal (a set of 
MOLS) if the squares in the set are pairwise orthogonal. Latin squares have been 
generalized to allow repetitions of elements in each row and column. 

Definition 3.1 Let n = m\. An F{n;X) frequency square is an n x n array in 
which each of m distinct symbols occurs exactly A times in each row and column. 
Moreover two such squares are orthogonal if when superimposed, each of the m? 
possible ordered pairs occurs times. 

The following result in fact contains Proposition 1.2 of 4j as a special case. 

Theorem 3.2 If there are E mutually orthogonal frequency squares of type F{n\ A) 
where n = mX, then M\{nX,nX — A^) > rriE. In particular, if q is a prime power 
and i > 1 is a positive integer, then 

M,.-.{q''-\ q'^'' - q''-') > q{q' - if/iq - 1). 

Further if i = 1, Mi{q, q — 1) = q{q — 1). 

Proof Label the rows and columns of each n x n frequency square by the elements 
0,1,...,?7, — 1. Then from each of the frequency squares, form a set of nA-tuples 
as follows. For each symbol i = 0, 1, ... ,m — 1, form an nA-tuple by listing the 
cell locations {k, I) where i occurs in the given square, proceeding row-by-row as k 
runs from to n — 1. Viewed as m blocks, each of size nX, of an affine resolvable 
design, these form a parallel class of size m. In total from the E squares, Em such 
nA-tuples are obtained, corresponding to E parallel classes. The entries of each nA- 
tuple are ordered pairs; form new nA-tuples by disregarding the first coordinate of 
each ordered pair. The resulting nA-tuples form the rows of an FPAx{nX, nX — X^). 
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For, since each symbol occurs A times in each column of a frequency square, each 
row of the array comprises A copies of each of the n column-headings. Any two rows 
of the FPA arising from the same parallel class will have distance n\. Any two rows 
derived from different classes will, due to the orthogonality of the corresponding 
frequency squares, agree in at most A^ positions, since agreement in p positions 
implies that some ordered pair occurs p times when the MOPS are juxtaposed. 
Hence the array has minimum distance nA — A^. □ 

For n = mX, it is known that the maximum number of mutually orthogonal 
frequency squares (MOPS) of the form F(n; A) is bounded above by (n — l)^/(m — 1). 
Further, if q is any prime power and z > 1 is a positive integer, then using linear 
polynomials in 2i variables over the finite field Fg, a complete set of F{q^; MOPS 
can be constructed. Specifically, take the polynomials aiXi+- ■ ■+a2iX2i where neither 
(ai, . . . , Oj) nor (aj+i, . . . , a2i) is the zero vector (0, . . . , 0) and no two of the vectors 
are nonzero Fg multiples of each other, i.e. (oi', . . . , 02/) 7^ e(ai, . . . , a2i) for any 
nonzero e G Fg. Further details may be found in Chapter 4 of [T^ . 

We remark in passing that, while the array obtained from Theorem 13. 21 is optimal 
in size when i = 1, it is not necessarily optimal for i > 1. This is in some sense 
expected because, in using these complete sets of mutually orthogonal frequency 
squares to construct error-correcting codes, the resulting codes are maximal distance 
separable only in the case when i = 1; see PU]- For example in the case q = i = 2, 
Theorem 14.61 vields an FPA2{8,4:) with more than 18 rows (see Example 14. 7|) . 

Another way to build frequency permutation arrays utilises finite fields, and may 
be considered as extending the approach of Theorem 2.4 of |3]. 

Theorem 3.3 Let L{x) = J^ll^o'^sX'^' ^ -^q^H- Denote by the degree of L{x), 
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and by r the rank of the matrix 



ML) 



( 



ai a, 
0L2 ot\ a 




,1-1 



,1-1 



V ^1-2 "Ls 




SO that 1 < r < i. Let < d < q 



Then 



d 



M,.-r{q\q^-dq')>Y, 



ii—r 



where Nj{q^) denotes the number of permutation polynomials over F^i of degree j . 



Proof It is a well-known result (see p 361 of ^S]) that the linearized polynomial 
L{x) is a permutation polynomial of F^i if and only if the determinant of the matrix 
A{L) is non-zero. More generally, the value set of L has cardinality q^, where r is 
the rank of A{L). So the linear transformation on Fgi defined by the polynomial 
L{x) has image of cardinality q^ and kernel of cardinality q^~^. Note that g*"** < qK 
Form an array as follows: for each permutation polynomial f{x) over Fg,, form 
a row by taking the images of the function L{f{x)) as x runs through the elements 
of the field Fgt. Each row is a A-permutation of length on m = q"^ symbols, each 
occurring with frequency A = q^~^. If f{x) and g{x) are permutation polynomials 
over Fgi of degrees at most d, then the polynomial L{f{x)) — L{g{x)) has degree at 
most dqK Hence (unless it is the zero polynomial) it has at most dq'' roots in Fgi, 
and so appropriately chosen f{x) and g{x) yield distinct rows of distance at least 
q^ — dqK We must now ensure that L{f — g) is not the zero polynomial. This happens 
if and only if the value set of the polynomial f — g lies wholly within the kernel of L, 
which has cardinality q^~^. Suppose first that f — g is non-constant. Now, f — g has 
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degree at most d < and, since a polynomial of degree d cannot have more than 
d roots in a field, its value set has cardinality at least [^-j^J + 1 > > g*"*". So the 
value set oi f — g cannot be contained entirely within the set of g*"*^ values mapped 
by L to zero, and hence L{f — g) is not the zero polynomial. For the constant case 
note that, for any permutation polynomial f{x), all f{x) + c with c G Fgi are also 
permutation polynomials. For f{x) + c to yield distinct rows, c must run through 
precisely one representative for each coset of the kernel of L; there are of these. 
Taking |j of the total number of permutation polynomials yields the desired number 
of rows. □ 

Observe that, in Theorem 13.31 if we take L to be the permutation polynomial 
x'^' \ we have maximal rank r = i and degree = q^^^, so we obtain a PA{q^, — 
dq^~^) of size ^^=1 A^j(g*)- 

To build an FPA with desired parameters, appropriate linearized polynomials 
may be chosen, whose properties are known in advance. The following corollaries 
illustrate two such constructions. Recall that the trace function TR : Fqi Fq is 
defined for a G Fgi by TR{a) = a + a'^ + a'^^ + ■ ■ ■ + a'^* \ More generally, letting 
i = gh and setting E = Fqi and F = Fgh, the trace function TRe/f '■ E F is 
defined for a G i^^ by 

TRE/F{a) + + a^'" + ... + a"?''"'". 

Corollary 3.4 Let q he a prime power and let i and h he positive integers such that 
h divides i. Let < d < q^. Then 

M,.M<l\l'-dq"')>J2iBr- 

Proof Let i = gh, let E = Fgi and F = F^h . Take L to be the generalized trace 
function TRe/f defined above; its kernel has cardinality q^'^ and its value set has 
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cardinality g'*. For any two permutation polynomials /, g over F^i of degree at most 

d < cl\ the value set of (non-constant) f — g has cardinality at least [^-^J +1 > (f ^^ 
and so TR{f — g) is not the zero polynomial. As in the proof of Theorem 13 . 31 dividing 
by q"^^^ deals with the case when f — g is constant. □ 
In the next section we consider how a PA may be converted into an FPA by 
appropriate substitutions on its symbols. If g is a prime power, a natural choice 
might be to apply the trace function to the rows of a PA{q^,d). However if, for 
example, there are two rows in the PA{q^, d) which differ by a constant a G Fqi with 
TR{a) = 0, then the resulting two rows in the FPA will be identical and so the 
rows will have distance 0. Thus applying the trace function to the elements of an 
arbitrary PA does not appear to be a good method to apply in a general setting. 

Corollary 3.5 Let q be a prime power and let i and n be positive integers such that 
n{< i) divides i. Let < d < g*~". Then 

Mq.{q\q^-dq-)>J2^^^ 

where Nj{q^) denotes the number of permutation polynomials over Fqi of degree j . 

Proof Let L{x) = x'^" — x in Theorem 13. 31 its roots are precisely the elements of Fqn. 
The polynomial L defines a linear transformation on Fqi whose kernel is the subfield 
Fqn and whose value set has cardinality g*~". For permutation polynomials f,g of 
degree at most d < g'""", non-constant f — g has value set of cardinality at least 
[^-^J + 1 > g" and so identical rows can arise only in the case when g(x) = f{x) + c 
with c G Fqn. □ 
We refer to jH] for a method for computing the value of Nj (g*) for any prime power 
g and positive integers i and j. Note, however, that the result of [S] requires consider- 
able computation to compute, and that the corresponding permutation polynomials 
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which arise from solutions to the system of equations in jH] must be constructed 
before the FPA can be built. 

We next indicate how affine resolvable designs can be used to construct FPAs. 
A balanced incomplete block design consists of a finite set V of v points, and a 
collection B of equally sized subsets of V called blocks, each of size k, such that 
every pair of distinct points of V occurs in exactly A blocks. A resolvable design 
has the additional property that the collection B of blocks can be partitioned into 
parallel classes (or resolution classes), such that every point of V occurs exactly 
once in each parallel class. An affine resolvable design{AIiD) is a resolvable design 
with the further property that any two non-parallel blocks intersect in precisely n 
points, where /i = ^ G N. When jj. = 1, the ARD is an affine plane of order k'^. 
Given an ARD, by labelling the blocks in each class and listing the blocks in which 
each point lies, an FPA may be constructed. 

Theorem 3.6 (i) Given an affine resolvable {v, k, A) design with r parallel classes, 
an FPAk{v, v — k) may be constructed of size r. 

(a) If there exist m MOLS of order n, then an FPAn{n^, n^—n) may be constructed 
of size m + 2. In particular, if q is a prime power, an FPAq^q^ , — q) may 
be constructed of size q+ 1. 

More details of this approach, including a proof of Theorem l3.6f i). may be found 
in [Sj. In the MOLS standard construction may be used to build m + 2 

parallel classes of an affine plane from the m MOLS, and these classes then used in 
part (i) to form an F P Anln"^ , n"^ — n). Equivalently, this FPA may be constructed 
directly by writing the rows of each of the m + 2 latin squares side-by-side to form 
new rows of length n^; it is clear from latin square properties that the resulting 
array is an FPAn{n^i — n). 
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Example 3.7 Using the following 2 MOLS of order 3 



L 



1 — 



1 2 

1 2 ,L2 

2 1 



1 2 
2 1 
12 



yields the following FPy43(9, 6); 



1 1 1 2 2 2 

1 2 1 2 1 2 

01212020 1 

1 2 2 1 1 2 0. 



Definition 3.8 An orthogonal array of size v, with r constraints, s levels and 
strength t, denoted OA[v, r, s,t\, is an r ^ v array with entries from a set of s > 2 
symbols, having the property that in every txv suhmatrix, every t x 1 column vector 
appears the same number ^ of times. 

The frequency array constructed in Theorem I3.(il is in fact an orthogonal array 
of strength 2. This gives rise to the following observation. 

Proposition 3.9 Every orthogonal array OA[v,r, s, 2] of strength 2 is an FPA2i{v,v— 
^) of size r. 

Proof In any row, each of the s symbols occurs with frequency ^. For any pair of 
rows, each of the pairs {i,j) of elements occurs times. In particular, each of 
the s pairs (i, i) occurs ^ times, and hence two rows agree pairwise in precisely ^ 
positions. □ 

Note that the FPAs obtained in this way are equidistant, in the sense that any 
two rows have distance precisely v ~ ^. For constructions of orthogonal arrays, see 
for example [6j; their connection with affine resolvable designs is explored in JT|. 

We end the section with a construction of frequency arrays from MDS codes. 
Recall that a g-ary (n, k) code is said to be maximal distance separable (MDS) if it 
satisfies the Singleton bound with equality, i.e. if d = n — k + 1. 
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Theorem 3.10 Given an [n, k, d] MDS linear code C over Fg, the array formed by 
taking the codewords of C as columns is an FPAqk-i{q^,q^~^{q — 1)). 

Proof Let C be an [n, fc, d] MDS linear code over Fq. Let G be a x n generator 
matrix for the code C, and write G = [G1G2 ■ ■ ■ Gn], where the Gi are the columns 
of G. Form an n x array A by taking the codewords of G as the columns of A. 
These are given by G^x^ = G^ as x runs through F^; the rows Gi^, . . . , Gn^ of G^ 
can be viewed as generating the rows of A. 

Each element of Fg occurs in each row of A with frequency q''~^, i.e. occurs g'^"^ 
times as the ith coordinate of the codewords of G. Let gi,...,gk be the elements 
in the i-th column Gi of G, and consider the equation aigi + ... + a^gk = b, where 
b,ai, . . . ,ak G Fg. Since Gi has at least one nonzero value, say in the j-th row, we 
can isolate the term ajgj = b — '^i^jO'igi- Then we can arbitrarily assign q values 
to each of k — 1 remaining a's, and uniquely solve the equation for aj since gj 7^ 0. 
Thus there are q^~^ solutions for each value of b in the i-th coordinate. 

Consider the distance between the two rows of the FPA corresponding to Gi^ and 
Gj^. We have the system of equations Gi'^-{xi, . . . , Xk) = a and Gj^-{xi, . . . , Xk) = (3 
{a, (3 G Fg). For an MDS code, any k columns (in particular, any two columns) 
of the generator matrix are linearly independent. Since Gi^ and Gj^ are linearly 
independent, this system of two linear equations in k variables will have rank 2, and 
thus solutions. This means that every ordered pair {a, (3) occurs q^'"^ times. 
Thus, in particular, the q ordered pairs {a, a), a G Fg are obtained g'^"^ times, so A 



Observe that this provides a direct construction for a class of FPAs described in 
Proposition 13.91 with v = q^ and s = q. Example 13. 71 mav alternatively be obtained 
by the MDS construction using the generator matrix 



is an FPA with distance q^ — qq^ = q^ — q' 



,fc-i 



□ 




1 

1 1 
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4 Constructing new FPAs from old 

In this section, we explore how one or more FPAs may be used as ingredients in the 
construction of new FPAs. 

Theorem 4.1 (i) Given an FPAx{n, d) of size N, a PA{n, d) may he constructed 
of size XN. In particular, Mx{n,d) < 

(a) Let I divide A. Given an FPAx{n,d) of size N, an FPAi{n,d) may be con- 
structed, of size jN . In particular, Mx{n,d) < jMi{n,d). 

Proof (i) Denote the FPAx{n, d) by A] let the symbol set of A be {0, 1, . . . , m — 1}. 
Using appropriate substitutions, A can be converted to a PA{n,d)^ A', of size A^. 
For a row R of A, moving from left to right, replace the A occurrences of a given 
symbol s by the sequence sA + 1, sA + 2, . . . , (s + 1)A (0 < s < m — 1). The new 
row R' is a permutation of 1, 2, . . . , ra. Since agreement between any two rows of 
A' can occur only at positions of agreement between the corresponding rows of A, 
the PA A' has minimal distance d. Now perform a cyclic shift on the entries of 
each substitution set {sA + 1, sA + 2, . . . , (s + 1)A} (0 < s < m — 1). This process 
can be repeated A times, to obtain A different substitutions for i?; all have pairwise 
distance n. Apply this process to each row of A\ the distance between new rows 
corresponding to different rows of A is at least d. Hence we have a PA{n, d) of size 
AA^. 

(ii) The proof is analogous to that of part (i). In this case, the substitution set 
for a given symbol s of the FPAx{n,d) comprises / copies each of a symbols. The 
generalization of the A cyclic shifts applied to the substitution sets, is the set of j 
permutations comprising an FPAi{X, A), described in part (ii) of Theorem 12.51 □ 

Example 4.2 An -FPy46(12,6) is constructed in Example 15.,'^ By Theorem \4 1[ 
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this FPA can he converted first to an FPA2,{12,Q)) and then to a PA{12,Q). We 
illustrate the use of the substitutions (without the cyclic shifts) on four sample rows. 



The first four rows of the FPAq{12, 6) are 

101011100010 

100101110001 
110010111000 
101001011100 



After substitutions, four rows of the FPA3{12, 6) are 

313134412242 
311313442224 

331131444222 
313113244422 

After substitutions, four rows of the FPAq{12, 6) are 

71829 10 11 3 4 5 12 6 

712839 10 11 4 5 6 12 

781293 10 11 12 4 5 6 

718239 4 10 11 12 5 6 

Converting a PA to an FPA by substitution is less straightforward in generaL 
The next result applies, for example, to an FPA arising from an orthogonal array. 

Proposition 4.3 Let n — mX. Let A be an FPAx{n, d) such that, between any two 
rows, each of the rr? pairs occurs precisely t times. Then A may be converted, 
by reduction mod r (where r\m) to an FPA]i(n,n — 

Proof Reduce the entries of A mod r. Each row of the new array is a A-permutation 
on r symbols with frequency ^. For any two rows in the new FPA, the pair of entries 
(i mod r, J mod r) agree, for each of the y values of j in the congruence class of i. 
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This yields ^ pairs for a given value of i, yielding such pairs in total, i.e. a 
minimal distance of — □ 

r 

The substitution technique may also be used on permutation arrays which have 
been constructed from latin squares. For example, given a PA obtained from The- 
orem 13.21 using a set of — 1 MOLS of order q\ applying the (g* — l)/(g — 1) 
substitutions from Theorem 9.20 of J3] to its entries, yields an FPA as described in 
the second part of Theorem 13.21 This approach allows the FPA to be built without 
constructing the corresponding sets of MOFS. 

In the proof of Theorem 14. ![ a pairwise distance of n is imposed on the set of 
new rows derived from any given original row. Relaxing this condition to minimum 
distance d, the A-cycle (or its frequency analogue) may be replaced by an appropriate 
(frequency) permutation array. This observation underlies the next result. 

Theorem 4.4 Let n = mX, and let Fi, . . . , Fi, be b FPA\{n, d) 's (not necessarily 
different). Let C be an FPAn{bn, c) of size N , where c > bd. Then an FPA\{bn, bd) 
may be constructed, of size iVmini<j<ft|Fj| . 

Proof Relabelling if necessary, construct Fi, ... ,Fb on disjoint symbol sets, so there 
are bn symbols in total. For each row in C, use the entries of the row as column 
headings, and place the n columns of each F^ under the n occurrences of the symbol 
i. The resulting array is an FPAx{bn,bd), of size mini<j<b|Fj|. Take the union 
of the arrays arising from each row of C to obtain an FPA of size A^mini<j<;,|Fj|. 
Agreement between rows of this FPA corresponding to different rows of C can occur 
only at positions where the rows of C agree, since the symbol sets are disjoint. There 
are at most c such positions, so any two rows of the new FPA have distance at least 
c > bd, and the array is an FPAx{bn, bd). □ 
A useful tool in building new arrays from old is the direct product. 
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Proposition 4.5 Let Xi be an FPAx{ni,a) of size Ni and let X2 be an FPAx{n2,b) 
of size A^2- Then an FPAx{ni + ^2, min(a, b)) may be constructed of size NiN2. 

In particular, for even n, given two FPAx{n, |), of sizes Ni and N2 respectively, 
an FPAx{2n, |) may be constructed of size N1N2, so that 

Mx{2n,-)>Mx{n,-f. 

Proof Relabelling if necessary, construct Xi and X2 on disjoint symbol sets, giving 
"^j^"^ symbols in total. Take the direct product of Xi and X2, i.e. Y = {{u,v) : 
u G Xi,v G X2}, where an ordered pair of codewords is interpreted as their con- 
catenation. Now, y is a set of A-permutations of length rii + ^2, with frequency A. 
Any pair of A-permutations in Y differ in at least min(a, b) positions, hence Y is an 
FPAx{ni + n2,in^iTi{a,b)). □ 
In j7], a permutation array is defined to be r -separable if it is a disjoint union 
of r Py4(n, n)'s of size n. We constructed an example of such a PA in part (i) of 
Theorem 14.11 We use this notion of a separable PA, i.e. a PA which is a disjoint 
union of other PA's, in the next result. 

Theorem 4.6 (i) Given a separable PA{n,d) which is the disjoint union of r 
PA{n, S) 's, each of size N, where 2d > 6, an FPA2{2n, 6) of size vN"^ may be 
constructed. 

(a) Given r MOLS of order n, an FPA2{2n,n) of size rrt^ may be constructed. If 
n is a prime power, an FPA2{2n,n) of size {n — l)n^ is obtained. 

Proof Denote the r PA{n, 5)'s by Fi, . . . , F^. For each i = 1, . . . ,r, form the direct 
product of Fj with itself, i.e. Zi = {{u,v) : u,v & Fj}. Then Zi is a set of A^^ 
A-permutations of length 2n, on n symbols, with frequency A = 2, and minimum 
distance 6. Take the union Z = ZiU . . .U Z.^. The A-permutations from different Zj 
have pairwise distance 2d > 6, and hence Z is an FPA2{2n, 6) of size vN"^. 
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By a result established in [7] and reproved constructively in ^2]; f MOLS of 
order n may be used to construct an r-separable PA{n,n — 1). When used in the 
above construction, this yields an FPA with 6 = n and 2d = 2n — 2 {> n for n > 2), 
i.e. an FPA2{2n,n) of size rn^. The last part follows by noting that, for a prime 
power n, a complete set of n — 1 MOLS of order n is obtainable. □ 

Example 4.7 By the construction from part (ii) of Theorem \4.(^ an FPA2{8,4:) of 
size 48 may be obtained from 3 MOLS of order 4. For example, the MOLS 



0123 0123 0123 

1032 2301 3210 
2301'^~3210'^~1032 

3210 1032 2301 



yield an FPA whose first 8 rows are listed below. 

1 2 3 1 2 3 

1 2 3 1 3 2 

1 2 3 2 3 1 

1 2 3 3 2 1 

1 3 2 1 2 3 
1 3 2 1 3 2 
1 3 2 2 3 1 
1 3 2 3 2 1 

The next two results generalize the direct product construction of Theorem 14.61 
A similar approach is explored in [S], in the context of constant composition codes; 
the reader is referred to [S] for more details, and for proofs of Theorem 14.81 and 
Theorem 14.91 

Theorem 4.8 Fori = 1, . . . ,b, let Xi be a separable FPAx{n, di) which is a disjoint 
union of ri FPAx{n, 6i) 's, r*i, . . . , FV^, with 'Y^di> min{5j}. Denote min{5j} by 5, 
andm.m.{ri} byr. Then an FPAh\{bn, 5) may be constructed, of size ^j=i(ni=i 1^/1)- 
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The direct product construction in Theorem l4.6l and Theorem l4.8l mav be adapted 
in various ways by choosing some subset of the direct product which has special 
properties. In [5, a recursive construction of PA's is given, which uses transversal 
packings; the next result indicates one way in which transversal packings may be 
used to construct an FPA from separable PAs. This construction may be applied to 
a set of disjoint separable PA's such as those obtainable from the MOLS construction 
of PHI, or to a single such PA with its subarrays permuted appropriately. 

Theorem 4.9 Let Xi, . . . ,Xk be k separable PA's, such that each Xi = PA{n,di) 
is a disjoint union of r PA{n, d/) 's, Ti^^\ . . . , Ti^'^\ and the Fj^-^-' 's may be ordered 
such that Fi'--'^ . . . , Ffc*--''' are disjoint for each j. Suppose there exist transversal 
packings Ti, . . . ,Tr , where each Tj has distance 6 and type iFi^-'-'l . . . iFfc^-'-'l. Denote 
di + ■ ■ ■ + dk by D, and denote the smallest sum of any 6 of the d/ by t. Then an 
FPAk{kn,d) may be constructed, of size Yl]=i l^jL where d = mm{t,D). 

We note that Theorem 3.2 of [4 may also be generalized to construct an FPAx{n, d) 
from k separable FPAx{ni, di)'s {1 < i < k). Replacing the PA(ni,(ij)'s by the 
equivalent FPA's in the statement and proof of this result, an immediate general- 
ization for A > 1 is obtained. 

5 Special cases 

An FPAx{n,d), where n = mX, may be viewed as an m-ary code with constant 
weight composition (A, . . . , A). In certain special cases, known results for constant 
weight codes provide bounds and constructions of relevance to FPAs. 

Proposition 5.1 If n = 2X, then an FPAx{n,d) of size M is a binary code 
{n,M,d) of length n, minimum (Hammimg) distance d and constant weight A. 
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In constructions and bounds are given for A{n,d,w), the maximum possible 
number of binary vectors of length n, Hamming distance at least d and constant 
weight w, for values of n up to 28. Observe that A{n, d, |) = Mn (n, d). The exact 
value of A{n, d, w), and corresponding constructions, is known for all lengths n < 11. 
If d is odd, then M^{n, d) = M;i{n, d+1), so only even distances need be considered. 
The following FPAs may be directly constructed, by use of Hadamard matrices and 
Steiner systems ([2]). 

Recall that a Hadamard matrix is a square matrix with entries +1, —1 whose 
rows are mutually orthogonal. Hadamard matrices of order n can only exist for 
n = 1,2 and n = 4k; it is conjectured that they exist for each n = 4k. Hadamard 
matrix constructions and properties may be found in Section IV. 24 of [U]. 

Theorem 5.2 (Theorem 10, Mii{n, |) = 2n — 2 and only if a Hadamard 
matrix Hn of order n > 1 exists. 

An (optimal) FPAn(^n, ^) may be constructed from the Hadamard matrix Hn 
as follows. First convert the entries of the 'half-frame' of +l's bordering if„ into 
— I's. Now take the non-initial rows of if„ and —H^, and convert the entries +1 to 
and —1 to 1 in every row. 

Example 5.3 Using the Hadamard matrix of order 12 (unique up to isomorphism) 
gives an FPAq{12, 6) of size 22. We list the first few rows. 

101011100010 
100101110001 
110010111000 
101001011100 

Combining Theorem 15 . 21 with Proposition 14 . 51 we see that, if n is an even number 
such that a Hadamard matrix of order n exists, then M^{2n, ^) > (2n — 2)^. For 
example, M6(24,6) > 484. 
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A Steiner system S{t, k, v) is a t — (f , k, 1) design, that is, a collection of fc-subsets 
(called blocks) of a f-set such that each t-tuple of elements of this f-set is contained 
in a unique block. When t = 3 and = 4, this called a Steiner quadruple system. 

Example 5.4 Using the Steiner quadruple system S'(3,4,8), an FPA4(8,4) of size 
14 may he constructed. The extended cyclic code {(1011000)1,(0100111)0} is one 
example; the code is constructed by taking cyclic developments of the vectors in 
parenthesis. 

We conclude by remarking that, in the study of PA{n, d) arrays, one builds the 
rows of the array by using permutations on n symbols, and in FPAx{n., d) arrays, one 
builds rows by using m distinct symbols, each repeated exactly A times. However, 
there is in fact no need for such uniformity of frequency, and one could consider the 
following, very general, setting. 

Let = Ai + ■ ■ ■ + Ar- be a partition of n. Then one could consider constructing 
arrays with the property that in each row, for i = 1, . . . ,r, the symbol i occurs 
exactly Aj times. From papers such as jl], there is motivation for studying such 
a general setting; in fact the corresponding constant composition codes have been 
widely studied; see for example |2]. Sets of F{n] Ai, . . . , A^) orthogonal frequency 
squares have been studied (see Chapter 4 in PJ'j)- However, we do not consider 
frequency permutation arrays with an arbitrary frequency vector n = Ai + ■ ■ ■ + A^ 
in this paper. 
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